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PACS 47.63. Gd “ swimming of microorganisms 
PACS 87.18. Gh - cell-cell communication 
PACS 87.17.Jj - Cell locomotion,chemotaxis 

Abstract -Bacteria suspension exhibits a wide range of collective phenomena arsing from interac¬ 
tions between individual cells. Here we show Serratia marcescens cells near an air-liquid interface 
spontaneously aggregate into dynamic clusters through surface-mediated hydrodynamic interac¬ 
tions. These long-lived clusters translate randomly and rotate in the counter-clockwise direction; 
they continuously evolve, merge with others and split into smaller ones. Measurements indicate 
that long-ranged hydrodynamic interactions have strong influences on cluster properties. Bacterial 
clusters change material and fluid transport near the interface and hence may have environmental 
and biological consequences. 
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Active systems are composed of self-propelled particles 
that can produce motion by taking in and dissipating en¬ 
ergy m- Examples exist at different length scales, from 
bacteria suspension [SHin] to flocks of birds [IIHIS]- Be¬ 
ing far from thermal equilibrium, active systems are not 
subject to thermodynamic constraints, such as detailed 
balance or fluctuation-dissipation theorem [HHTB] . This 
renders the physics of active systems much richer than that 
of thermal systems. For example, collective motion with 
extended spatio-temporal coherence has been reported in 
many active systems [MsjinHiiinHis]- Such coherent 
motion can arise from local interactions that align a par¬ 
ticle’s motion with its neighbors through biological coor¬ 
dination [nma or physical interactions [mils]- 

Active systems without alignment interactions also ex¬ 
hibit interesting collective behavior. Theoretical models 
have shown that systems with a density-dependent motil¬ 
ity phase separate into dense dynamic clusters and a di¬ 
lute gas phase Numerical simulations of repul¬ 

sive self-propelled disks confirmed the theoretical predic¬ 
tion of phase separation [2T1[^. Effects of motility, at¬ 
tractive interaction, and hydrodynamic forces have been 
extensively explored in simulations |23fBB] . On the ex¬ 
perimental side, dynamic clusters have been observed in 
Janus particles (platinum-coated [57] and Carbon-coated 
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[55] ') and colloidal particles with an embedded hematic 
cube [55]. Schwarz-Linek et al observed clusters of motile 
bacteria when they added polymers to bacteria suspension 
to induce depletion attraction between bacteria [SO]. In a 
very recent paper m, Petroff et al reported that Thiovu- 
lum majus bacteria form two-dimensional crystals near a 
liquid-solid interface. Understanding the origins and prop¬ 
erties of these dynamic clusters may provide new insights 
into emergent behaviors of active matters and open up 
possibilities to build novel materials [15] . 

In this letter, we report experimental results for a new 
type of bacterial clusters formed near an air-liquid inter¬ 
face in a pure suspension without depletant agents. Fluid 
dynamic calculation and flow visualization are used to 
show surface-mediated hydrodynamic interactions can ex¬ 
plain the formation of these clusters. We further quantify 
the statistical and dynamic properties of bacterial clusters 
and show long-ranged hydrodynamic forces have impor¬ 
tant influences on cluster properties. We conclude with 
discussions on related research and on possible technolog¬ 
ical and environmental implications of our work. 

Experiments - Our experiments are carried out in 
drops of wild-type Serratia marcescens (ATCC 274) bacte¬ 
ria, which are propelled by a bundle of a few rotating flag¬ 
ella m- For cultivation, small amount of bacteria from 
frozen stock is put in 4 ml of Luria Broth (LB) growth 
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Fig. 1: (color online), (a) Dynamic clusters observed in a sam¬ 
ple with bacterial density 0 = 0.058^m“^. (b) Clusters dis¬ 
solve after flagellar motor is damaged by strong light. The 
false colors in (a-b) show the time evolution of bacteria origi¬ 
nally belonging to different clusters; image in (b) is taken 0.3 s 
after that in (a), (c-d) show a merging event of two clusters in 
a sample with cj) = 0.044/rm“^. Overlayed arrows show cell ori¬ 
entations. Two ellipses are drawn to highlight difference in cell 
motility in/outside clusters, (e) Probability density functions 
of averaged angular velocity of bacteria in(red) /outside(green) 
clusters, (f) Schematic of the experimental setup (not to scale). 
Scale bars in (a-b) and (c-d) correspond to 20 /rm and 15 ^tm, 
respectively. A coordinate frame is defined in (a) and (f): the 
z axis points into the bulk fluid from the trapping (xy) plane. 


medium consisting of 0.5% yeast extract (Sangon G0961), 
1% Tryptone (Sangon TN5250), and 1% NaCl (Sigma). 
Bacteria is incubated for 13 hours to a stationary phase in 
a shaking incubator which operates at 30° C and 200 rpm 
shaking speed. We then extract 1 ml bacteria solution and 
re-grow bacteria in 10 ml fresh LB growth medium supple¬ 
mented with 5^g/ml of A22 for 4 hours at 33°C and 200 
rpm. A22 is a small molecular inhibitor of a protein MreB, 
which is needed to maintain rod-like shape of many bacte¬ 
ria m- s. marcescens cells growing in a media with A22 
are short ellipses with a mean aspect ratio of 1.2. The fi¬ 
nal bacteria solution is further diluted in distilled water to 
generate samples with various bacteria densities. To con¬ 
trol bacterial motility, l//g/ml photosensitizer FM 4-64 is 
added to the bacterial suspension m- In the presence of 
FM 4-64, the bacteria will be temporally paralyzed upon 



Fig. 2: (color online). Computed flow field around a model 
bacterium on two planes: y = 0 in (a) and z = 0.54 fim in 
(b). Color represents magnitude of velocity projection in the 
plane and arrows denote the direction of flow. Profiles of ra¬ 
dial (blue) and tangential (red) velocity components along the 
white dashed line in (b) are plotted as symbols in the insert 
with fits (lines) to analytical expressions derived in Supple¬ 
mentary Text ST.pdf. Computation is carried out with the 
Regularized Stockslet method |35II36| and details can be found 
in Supplementary Text ST.pdf. 


exposure to strong light. A lOOW Mercury lamp(Nikon 
C-SHGl) is used to activate photodynamic effects. 

Sample is enclosed in a sealed chamber (cf. Fig. 1(f)) 
which consists of a silicone well, a plastic spacer, and 
two cover slips (0.13 mm thick). The well has a diam¬ 
eter of 0.8 cm and a height of 0.1 cm. In some experi¬ 
ments, 1 ^m super-paramagnetic tracer beads are added 
for flow visualization. Beads are confined to the inter¬ 
face by a permanent magnet and the strength of confine¬ 
ment can be tuned by moving the magnet relative to the 
sample. Two-dimensional bacterial or tracer motion in 
the trapping plane is imaged through a 60X phase con¬ 
trast objective (Nikon ELWD ADL 60XC) and recorded 
by a camera (Basler acA2040-180km). We use a holo¬ 
graphic microscope to measure three-dimensional motion 
of tracer particles and bacteria. A red LED is used for 
illumination and recorded holograms are analyzed by the 
Rayleigh-Sommerfeld back-propagation method to extract 
the spatial coordinates of the scatter [ 23135 ] ■ 

Bacteria clusters - S. marcescens is known to ad¬ 
here strongly to the air-liquid interface possibly due to 
hydrophobic interactions [SM]- We use a holographic 
microscope to record such attaching events. Supplemen- 
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Fig. 3: (color online), (a-c) Interaction between a tracer bead (blue +) and a bacterium in the trapping plane. The red arrows 
mark the instantaneous orientation of cell body, (d) Time series of angular velocity (blue) of the cell and the separation (green) 
between the cell and bead, (e-f) Three-dimensional flow around a bacterium (red *) visualized by a tracer bead (black sphere). 
In each panel (e-g), the black-and-white image in xy plane is the raw hologram and the blue line is the tracer trajectory, (h) 
Time series of the z-coordinate of the tracer (red) and the separation (green) between the bead and the cell. Arrows in (d) and 
(h) mark when data in other panels are recorded. 


tary movie SI.mp4 0 shows two typical events: bacteria 
swim from the bulk towards the interface and get trapped 
near the interface. Trapped bacteria can move freely in 
the trapping (xy) plane; their centers-of-mass are tracked 
to quantify translational motion. Rotational motion is 
quantified by following the principal axes of the elliptical 
bacteria. The main control parameter in experiments is 
bacteria density in the trapping plane which is quantified 
by the number of bacteria per unit area and denoted as (j>. 

As shown in Fig. l(a-b), trapped bacteria aggregate 
to form dynamic clusters which are identified through 
Voronoi analysis of bacterial positions (see Supplementary 
Text ST.pdf □ for details). Clustering formation is a ro¬ 
bust phenomenon and occurs under a wide range of bacte¬ 
ria densities and in various liquid environments including 
LB media, motility buffer, and water. However, formation 
of dynamic clusters requires bacterial motility. As shown 
in S2.mp4[l, clusters dissolve immediately after we reduce 
the motility by using strong light [31]. Clusters reappear 
when bacteria motility is partially recovered after the light 
irradiation is stopped. Heavy metal ions (CUSO 4 ) are also 
used to reduce motility [42]; similar results are obtained. 

To understand the connection between motility and 
cluster formation, we zoom-in and investigate how bac¬ 
teria move in a cluster. We mark a bacterium in a clus- 

^The movie SI.mp4 shows the bacterium trajectory as a blue 
line and the raw hologram as a black-and-white image. A black dot 
marks the instantaneous position of the bacterium. The objective 
is focused on the interface. All submitted movies are compressed in 
Xvid Codec and movies in MJPG format are available at: 
http://ins.sjtu.edu.cn/people/hpzhang/EPL/Movies.rar 

^ST.pdf contains detailed information of data analysis procedures 
and of fluid dynamic analysis and is available at: 
http://ins.sjtu.edu.cn/people/hpzhang/EPL/ST.pdf 

^The movie S2.mp4 shows identified clusters in colors on phase- 
contrast images. Light irradiation starts at 16.8 s in the video. 


ter with a red ellipse in Fig. l(c-d) (cf. S3.mp4 0 ); the 
marked bacterium rotates its body in the trapping plane 
with an angular velocity of = —20 rad/s. In contrast, 
a bacterium outside cluster (marked in green) shows lit¬ 
tle change in its body orientation. In Fig. 1(e) we plot 
the probability distribution functions of angular velocity, 
P{nz), for bacteria in and outside clusters. While P{ftz) 
peaks around zero for bacteria outside clusters, P{i}z) 
shows a second peak at -23 rad/s for bacteria in clusters. 
We next use fluid dynamic calculation and flow visual¬ 
ization to show that the bacteria with large can form 
clusters through hydrodynamic interactions. 

Hydrodynamic interactions - S. marcescens bacte¬ 
ria swim by rotating their flagellar bundles. When viewed 
from the front (body) of a bacterium, the bundle rotates 
clockwisely and the cell body rotates in the opposite direc¬ 
tion to achieve hydrodynamic torque balance. Fast body 
rotation in Fig. 1(e) suggests the bacteria in clusters ori¬ 
ent their bundles perpendicular to the interfac(0. Bacteria 
in such a configuration can generate fluid flow that leads 
to cluster formation |3TJ|43l|44] . To illustrate the mecha¬ 
nism, we numerically compute fluid flow around a bacterial 
model that is oriented perpendicular to the interface. As 
shown in Fig. 2(a), the model has a 1-^m-diameter spher¬ 
ical body whose center is located at (0, 0, 0.54 ym) and is 
driven by a rotating flagellum. The translational degrees 
of freedom are frozen for the model. The model bacterium 
exerts a force in the positive z direction on the fluid, draws 

^In the movie S3.mp4, phase-contrast images are shown in the 
background, cell orientations are marked by arrows, and ellipses are 
used to highlight two bacteria. 

®A hydrodynamic mechanism was proposed to explain perpendic¬ 
ular orientation of T. majus cells m- Such a mechanism only works 
at a no-slip boundary and doesn’t apply in the case of S. marcescens 
near a liquid-air interface. 
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Fig. 4: (color online). Statistic and dynamic properties of clus¬ 
ters measured for six different bacterial densities (color-coded 
according to the legend in (a)). Quantities measured include 
probability distribution function (a), mean translation speed 
(b) and mean rotation period (c). Solid lines in (b) and (c) are 
salings derived from force and torque balance (see text). 


body shows little rotation perpendicular to the interface, 
i.e. riz is small [51151H5]. 

Important flow features in Fig. 2 are qualitatively con¬ 
firmed in experiments with two types of tracer beads. The 
first kind is strongly confined in the z-direction and only 
probes fluid flow in the trapping plane. A typical result is 
shown in Fig. 3(a-d) and Movie S4.mp4[^. A tracer bead 
is drawn to the bacterium from t = 0tof = 0.5s, which 
demonstrates the inward radial flow. While the bead and 
bacterium are bound, they rotate around each other coun- 
terclockwisely, which is a manifestation of the tangential 
flow in Fig. 2(b). At t = 7 s, the bacterium stops rotating 
and swims away from the bead. In the second experi¬ 
ment, a tracer bead is subject to weak confinement and 
can be advected also in the z-direction. As shown in Fig. 
2(e-h) and Movie S5.mp4 0, the bead is drawn towards 
the bacterium in the trapping plane until t = 3.5 s, when 
the separation between the bead and bacterium is about 
1 /rm. The bead is then advected quickly into the bulk, 
which demonstrates strong flow into the bulk around the 
bacterium. 


Tracer particles are also used to visualize flow around 
bacteria clusters. Movie S6.mp4@ shows that clusters at¬ 
tract tracers and advect them into the bulk, mimicking 
results in Figs. 2 and 3. This supports the following pic¬ 
ture: bacteria in clusters orient their flagella perpendicu¬ 
lar to the interface; they generate inward radial flow that 
attracts neighbors to form clusters and counterclockwise 
tangential flow that drives clusters into rotation. 

Cluster properties - Bacteria clusters observed in ex¬ 
periments are highly dynamic; they constantly evolve and 
change their sizes. We defined cluster size as the num¬ 
ber of the constituent bacteria, n. Probability distribu¬ 
tion functions for finding a cluster of a given size at six 
bacteria densities </> are shown in Fig. 4(a) . As the den¬ 
sity increases, probability to find large clusters increases. 
Probability distribution function decays exponentially for 
large n. Similar exponential distributions have been ob¬ 
served in many previous studies [5115] and may be modeled 
by fusion-fission processes [46]. 


We compute the following quantities to quantify trans¬ 
lation and rotation of the /th cluster which contains 
the fth bacterium at a location fij and with a velocity 
Vij. The center-of-mass of the /th cluster is located at 


Ri = where {■)^ denotes an average over all n/ 

bacteria in the /th cluster. Speed of the center-of-mass 


is Vi 


Vi 


{vi,i)^ \ ■ The angular speed of the /th 


in fluid along the interface (at z = 0) and pushes fluid to 
the bulk. Rotation of the flagellum and cell body also pro¬ 
duces a tangential flow component that can be seen in Fig. 
2(b), especially in the region close to the cell body. In the 
far field, radial (W (r)) and tangential {Vg (r)) components 
decays as: W {r) oc r~^ and Vg (r) oc r~^, as shown in the 
insert. Bacteria outside clusters in our experiments likely 
orient their bundle parallel to the interface and the cell 


®In the movie S4.mp4, a passive tracer confined to the interface 
visualizes attractive flow around a rotating bacterium. 

^The Movie S5.mp4 shows three-dimensional motion of a passive 
tracer around a rotating bacterium. 

® Movie 6 shows three-dimensional motion of a passive tracer 
around bacterial clusters. A l/zm tracer bead (marked by red-h) 
is drawn towards a bacterial cluster during the first 1.7s, and then 
advected quickly into the bulk, as shown by the enlarging interfer¬ 
ence rings in the hologram. 
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cluster is defined as: w/ = 




whose 


inversion is the rotation period: Tj = Averaging Vi 
and Tf over all clusters of size n, we have mean trans¬ 
lation speed V (n) = (h/)^ and mean rotation period 

Results in Fig. 4 (b-c) show that, for a given (j), larger 
clusters translate and rotate slower than smaller ones. 
This dependence is qualitatively consistent with a simple 
model of active clusters made of self-propelled particles 
|31] . Each particle is driven by a propulsive force F^p in the 
Ci direction and has a drag coefficient of 7 ; they interact 
through pair-wise forces in the radial and tangential direc¬ 
tions . If n particles in a cluster are randomly oriented, to¬ 
tal propulsive force on the cluster is FspCi] = Fsp'nP'^. 
The total friction coefficient is the sum of all particles: 
ny. The average velocity of the center-of-mass is the ra¬ 
tio of the total force to the total friction coefficient and 
scales as: V oc n~^'^. Using a similar torque-balance argu¬ 
ment we can get a scaling law for rotation period: 

T (X n . These two scalings are shown in Fig. (b) and (c) 
as black lines; discrepancies with experimental results can 
be clearly seen. 

First, we notice that cluster motion in experiments de¬ 
pends strongly on global densities, (j). Fig. 4 (b-c) show 
that a cluster of a given size translates faster and rotates 
slower in a system with a higher density. The dependence 
on system density likely arises from long-ranged hydro- 
dynamic interactions which enable fluid disturbances to 
propagate far and couple different clusters over large dis¬ 
tances. Second, Fig. 4(b) shows that cluster velocity V (n) 
scales as U oc rather than V oc This is also 

possibly related to hydrodynamic interactions which de¬ 
termine how cluster friction and net propulsive force scale 
with cluster size. In the limit of Stokes drag, cluster fric¬ 
tion scales as y/n |48j in two dimensions, which suggests 
that net propulsive force of clusters may scale as 
Data in Fig. 4(b-c) indicate that hydrodynamic interac¬ 
tions in our experiments are too complex to be represented 
by simplified pair-wise forces,. 

Discussion - Inward radial flow in Fig. 2(a) is similar 
to feeding flow that many micro-organisms use to gather 
food from the fluid environment near an interface mM- 
The same flow pattern is also used to explain the forma¬ 
tion of bound Volvox pair |49| and to explain the attrac¬ 
tive force between thermophoretic colloids [501152) . Such 
an attractive boundary flow should exist near any low- 
Reynolds number swimmer that is oriented perpendicular 
to a fluid or solid boundary and swims into the boundary; 
consequently, these oriented swimmers experience effec¬ 
tive attraction and can be hydrodynamically assembled 
into clusters. This provides a new mechanism, besides 
phoretic m, depletion m interactions and self-trapping 
effects m, to generate clusters of active particles. 

Petroff et al E] recently reported that T. majus bacte¬ 
ria form two-dimensional crystals near a liquid-solid inter¬ 
face through a mechanism similar to that in S. marcescens 


clusters. However, T. majus and S. marcescens systems 
are significantly different in at least two aspects. First, 
the strength of attractive interaction is different. Petroff 
et al used the product of propulsive force and swimmer 
size to estimate the energy scale (denoted as E) for attrac¬ 
tion between cells. T. majus have an average diameter of 
8.5 /rm and swim at a speed of 600 /xm/s; it was found 
that the attractive energy is much larger than thermal 
energy: E ~ lO^fc^T. Consequently, T. majus crystals 
are very stable and can contain up to a thousand cells. 
On the other hand, S. marcescens, like many other com¬ 
monly studied bacteria [ziiiHiisn] and artificial swimmers 
ESlIlS], are approximately ten times smaller in both size 
and speed than T. majus\ the attractive energy scale is 
1000 times smaller and is on the order of lOfcsT. There¬ 
fore, fluctuations play a more important role in the S. 
marcescens system and render S. marcescens clusters dy¬ 
namic and constantly evolving. Second, T. majus crystals 
and S. marcescens clusters form near different hydrody¬ 
namic boundaries. A no-slip boundary makes the hydro- 
dynamic attraction, /, between two T. majus cells decays 
rapidly as their separation, r, increases: / ~ r“^; for 
S. marcescens clusters near a free-slip boundary, we have 
/ ^ r~^. Consequently, while hydrodynamic interaction 
between T. majus cells are severely screened, Fig. 4 shows 
that long-ranged hydrodynamic interactions influence S. 
marcescens clusters properties. 

Conclusion - In summary, we have investigated dy¬ 
namic clusters of S. marcescens bacteria near an air-liquid 
interface. Fluid dynamic calculation and flow visualiza¬ 
tion suggest that the constituent bacteria of these clusters 
orient their flagella perpendicular to the interface. Bac¬ 
teria in such a configuration generate radial flow that at¬ 
tracts neighbors to form clusters and tangential flow that 
sets clusters into counter-clockwise rotation. We measured 
statistical properties of bacteria clusters and showed clus¬ 
ter properties are affected by long-ranged hydrodynamic 
interactions. S. marcescens clusters efficiently change ma¬ 
terial and fluid transport near the air-liquid interface; they 
may have environmental and biological consequences |5]. 
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